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Abstract 

The free energy density of the SU(3) gauge theory at temperatures 
r/Tc = 4/3,3/2 and 2 is calculated on lattices with temporal extent as 
small as Nt = 2, 3 and spatial extent Ncr = 4Nt using parametrized 
fixed point actions. Although cut-off effects are seen, they are hugely 
suppressed with respect to Wilson and Symanzik-improved actions and 
at Nt = 3 there is already a good agreement with the continuum limit as 
extrapolated from the results with the Wilson action at Nt = 6 and 8. 



PACS number(s): 05.70.Ce, ll.lO.Wx, ll.lO.Hi, 12.38.Gc, 12.38.Mh 



1 Introduction 



Understanding the high temperature phase of QCD is of great importance es- 
pecially in astronomy and cosmology. In such a phase quarks and gluons are 
expected to behave like a plasma of weakly interacting particles approximately 
described by an ideal gas. To find signatures of this plasma is the central goal 
of present and forthcoming experiments in heavy-ion collisions at CERN and 
BNL. 

Owing to the well-known infrared problem Q of finite temperature QCD, 
the perturbative expansions of thermodynamic quantities can be calculated only 
up to order and their convergence is very bad even at temperatures much 
larger than |2[ |^. Lattice numerical simulations are therefore an essential 
tool in investigating QCD at high temperatures, where non-perturbative effects 
seem to be important at least in the vicinity of the phase transition. 

In lattice numerical simulations, the temperature T and the volume V are 
determined by the lattice size N'^ x Nt, < Na through 

T = ^ , V = {N^af , (1) 

and (anti-)periodic boundary conditions have to be imposed in the time di- 
rection for (Fermi) Bose fields. The lattice spacing a depends on the coupling 
f3 = 2N/g'^; the asymptotic form of this dependence can be known analytically 
from perturbation theory. 

At high temperatures, high momentum modes give relevant contributions to 
the free energy density. In this regime, thermodynamic quantities are therefore 
strongly influenced by the finite lattice cut-off introduced by the discretiza- 
tion of the theory. This cut-off effect manifests itself in O(a^) (i.e. 0{1/N^) 
at fixed T) corrections which vanish in the continuum limit. How large A'^^ has 
to be in order to allow reliable extrapolations to the continuum depends on the 
specific regularization of the theory, i.e. on the lattice action. 

The Wilson plaquette action Q is the simplest action which can be used in 
numerical simulations of the pure gauge theory. It has however the disappointing 
drawback of introducing large cut-off dependences in physical observables: in 
zero temperature simulations, the continuum limit is reached at lattice spacings 
smaller than 0.1 fm, so physics at the hadronic scales can be simulated only on 
uncomfortably large lattices. One can have a feeling of the relevance of cut-off 
effects with the SU(N) Wilson action looking at the expressions for the energy 
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density e and the pressure p of the gas of gluons in the T oo Umit 
e 3p 



. - . - (N' - 1) — 



(2) 



The deviation from the Stefan-Boltzmann law (iV^ — l)7r^/15 is very prominent: 
for Nr = 4, where also 0(l/iV^) contributions are not negligible, the deviation 
from the continuum is about 50% Q. The relevance of cut-off effects in SU(3) 
thermodynamic quantities at temperatures T equal to a few times Tc is less 
dramatic [0, ||, but large enough to call for lattices up to 32'^ x 8 before safely 
extrapolating to the continuum. Since fluctuations in thermodynamic quantities 
go like l/y/V4 ~ N~'^, a constant accuracy requires a computational effort which 
increases at least like N^'^ (here, a fixed Na/Nr ratio has been assumed). 

Adopting the Symanzik program Q to improve actions reduces indeed lat- 
tice artifacts. Since the improvement procedure is perturbative, the benefit is 
expected only at large values of the coupling /3, i.e. at large temperatures T for 
fixed Nt. Indeed, in SU(N) lattice gauge theory in the T — > oo limit, the first 
correction to the Stefan-Boltzmann law can be made 0{\/N^), if the action is 
improved at the tree level with alx2ora2x2 plaquette, or even 0{1/N^), if 
it is improved with a 2 x 2 and a 3 x 3 plaquette Q . For Nr = 4 the deviation 
from the Stefan-Boltzmann law is 1.34% with the 1x2 improvement and 8.77% 
with the 2 X 2 one @. 

As a radical solution of the problem of cut-off effects in lattice simulations, 
P. Hasenfratz and F. Niedermayer have recently suggested to use perfect 
lattice actions, i.e. actions which are completely free of cut-off effects. The 
existence of such actions is well established in the Wilson's renormalization 
group (RG) theory |l^. The fixed point (FP) action and the renormalized 
trajectory of a RG transformation define a perfect quantum action. For an 
asymptotically free theory, the FP action itself reproduces the properties of 
the continuum classical action - it is the perfect classical action. A formal 
argument of perturbative RG implies that FP actions are also 1-loop (quantum) 
perfect jlj, ^ . FP actions arc therefore the first step in the ambitious program 
of building perfect lattice actions. 

In Ref . [0 a general method for the determination of FP actions for asymp- 
totically free theories was proposed, which was successfully applied to the 2-d 
0(3) non-linear cr-model on the lattice. A parametrization of the FP action 
suitable for numerical simulations was found showing no cut-off effects up to 
correlation lengths as small as f ^ 3 lattice spacings, i.e. well away from the 
continuum limit. Moreover, a perturbative calculation of the mass gap in 
a periodic box of size L, m{L), using the lattice FP action showed that cut-off 
effects are absent also on the 1-loop level. 
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A similar program has been carried out for SU(3) lattice gauge theory [|l3[ ^ 
starting from two different RG transformations (type I and type II in the follow- 
ing). A few-parameter approximation has also been suggested | p5| for both type 
I and type II FP actions in terms of the plaquette and the twisted perimeter- 
six loop. Preliminary scaling tests of the type I FP action on the torelon mass 
and the static qq potential, using the critical deconfining temperature to set 
the scale, have already shown a clear improvement with respect to the Wilson 
action [T^ . 

The aim of this paper is to perform a study of the cut-off effects of pa- 
rametrized FP actions in bulk thermodynamic quantities of the SU(3) gauge 
theory at finite temperature. Indeed, the contribution of high excitations to 
these quantities is important, making lattice determinations quite sensitive to 
the cut-off. FP actions are implemented on lattices as small as 8'^ x 2 and 12'^ x 3 
to calculate the free energy density at temperatures T/Tc ~ 4/3, 3/2, 2. 

The paper is organized as follows: 
in Section 2 the construction of a (parametrized) FP action is shortly recalled 
and some comments are presented about the expected cut-off effects in the free 
energy density; 

in Section 3 elements and techniques of lattice finite temperature SU(3) ther- 
modynamics relevant to the remaining part of the paper are briefly reviewed; 
in Section 4 Monte Carlo (MC) determinations of the free energy density with 
parametrized FP actions are compared with previous results with Wilson and 
Symanzik-improved actions. 



2 Parametrized FP actions and cut— off effects 
in the free energy density 

The partition function of an SU(N) gauge theory defined on a hypercubic lattice 
is 

Z = J DU e-^^<^''^ , (3) 

where PA{U) is some lattice regularization of the continuum action expressed 
in terms of products of link variables U^{n) = e*'^^''"^ eSU(N) along arbitrary 
closed loops. If (3, 02,03, .. . denote the couplings of (3A{U), the action can be 
represented by a point in the infinite dimensional space of the couplings (see 
Fig. |l|). 

Under repeated RG transformations the action moves in the space of cou- 
plings. A RG transformation with scale factor 2 can be defined in the following 
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C2, . . . 




Figure 1: Schematic flow diagram for asymptotically free gauge theories. 



way: 

^-I3'A'{V) ^ I j^jj ^-f3(A(U)+T{U,V)) . (^4) 



here V is the coarse link living on the lattice with spacing a' = 2 a, related to 
a local average of the original link variables U ; T{U, V) is the blocking kernel 
which defines this average, normalized in order to keep the partition function 
(and therefore the free energy) invariant under the transformation. There is 
a wide freedom in defining the RG average; this freedom can be exploited for 
optimization in numerical simulations. 

On the critical surface f3 — oo, Eq. (Q) leads to the saddle point equation]^: 
A'{V) =mm[A{U)+T{U,V)] ; (5) 

the FP point of this transformation is therefore 

A^P{V) ^ mini A^P{U)+T{U,V)] . (6) 

An important consequence of this equation is that if the configuration {V} 
satisfies the FP classical equations and is a local minimum of A^^(y) then 

^/3' = /3 — 0(1) owing to asymptotic freedom. 
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the fine configuration {U{V)} minimizing the r.h.s. of Eq. (||) satisfies the 
FP classical equations as well and the value of the action remains unchanged, 
A^P{V) = A^P{U{V)) Q As a consequence, the FP action has exact 
scale invariant instanton solutions with the same action as in the continuum 
theory. 

Although the FP lies on the critical surface /3 = oo, by the expression "FP 
action" we mean (3A^^{U) which corresponds to the line which leaves orthog- 
onally the critical surface, i.e. I3A^P{U) P'[A^P(U) + 0{l/(3^)]. For very 
large /3, the FP action stays close to the RT, i.e. it is a good approximation of 
the perfect action. This follows from a formal argument |l3| of perturbative 
RG implying that FP actions are perfect not only in the classical limit /3 — oo, 
but also at 1-loop order of perturbation theory. The range of /3 for which the 
FP action and the RT stay close together cannot be known a priori just by 
looking at the RG transformation. 

The FP equation (^) is highly non-trivial. On smooth V configurations, it 
can be expanded in powers of the vectors potential and studied analytically. At 
the lowest order in the expansion of the vector fields the FP action is quadratic 
and can be found analytically ||T^: it describes massless free gluons with an 
exact relativistic spectrum for all values of the momenta, i.e. not restricted to 
the first Brillouin zone. 



On configurations with large fluctuations, numerical methods have to be 
used to solve Eq. (H). For practical reasons, however, it is useful to find a 
general parametrization for the FP action, simple enough to be used in Monte 
Carlo simulations. In Ref. ||l3| two RG transformations (type I and type II) 
were defined depending both on two parameters which were optimized in order 
to make the quadratic part of the FP action the most short ranged. The type 
I FP action (which is used in the following) was parametrized with powers of 
the real part of the trace of the plaquette and the twisted perimeter-six loop, 
calling for (only) 8 parameters. 

The source of cut-off effects in determinations of the free energy density on 
the lattice using parametrized FP actions as regularization is three-fold: 

i. Although the blocking kernel T{U, V) is normalized in order to keep the 
partition function invariant under the RG transformation, a field-independent 
term in A'{V) appears in Eq. (jj) as a result of integration on the fine links 
U. The amount of this constant term can be singled out plugging the flat 
configuration of the coarse lattice (V^(rtB) = 1 for every site ub) into Eq. (^: 
it is easy to see that it is given by the sum of quantum fluctuations around the 
flat configuration of the fine lattice. One usually discards this term since it does 
not affect the dynamics of the coarse field. For Nr, (a few times) larger than 
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the "size of the action" this term is proportional to the lattice volume N^Nr 
and gives again a constant (independent of Nr and No-) contribution to the 
free energy density which can be discarded as well. However, for very small Nr 
(and/or N^) the field-independent part will depend non-trivially on Nr- This 
produces an artifact which decreases exponentially with Nr. In practice, it will 
show up only for Nr — 2. 

It should be stressed that all the above observations hold even in the classical 
limit. 

ii. The perfect action at finite f3 is represented by the points on the RT. At 
moderate values of (3, the FP action deviates from the RT: this deviation causes 
cut-off effects. 

iii. The parametrization of a FP action introduces additional cut-off ef- 
fects. This is especially true for the few-coupling parametrizations needed to 
perform numerical simulations. In most cases parametrizations are reliable only 
in the region of correlation lengths in which the numerical procedure to fix the 
parameters was performed. 

While cut-off effects referred in ii. are not fully under control, those referred 
in iii. can be considerably reduced if a "good" parametrization is found. Short 
range actions are easier to parametrize than extended actions. Experience with 
RG suggests that a good averaging procedure in the block transformation will 
decrease the extension of the corresponding FP action. In Ref. |]l6| a new 
RG transformation has been proposed with a smoother and more rotational 
symmetric smearing kernel T(U,V), with a larger number of free parameters 
used to optimize the resulting FP action. Preliminary scaling test of this type 
III FP action with the static qq potential have already shown an improvement 
of the performances with respect to the type I FP action (using the same 8- 
coupling parametrization) |l6| . 



3 SU(3) thermodynamics on the lattice 



All the continuum thermodynamic quantities for a system with zero chemical 
potential can be derived from the free energy density]^ 

/ - -^lnZ{T,V) , (7) 

where V is the volume of the system, T the temperature and Z the canonical 
partition function. The energy density e and the pressure p are determined by 

^The Boltzmann constant is set to one. 
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the formulae 

T^ d\nZ{T,V) 

' Y — Wf — ' 

. ^ T^l™ . (9) 

For homogeneous, large {N^ > i) systems 91nZ(T, V)/dV = \nZ{T, V)/V, so 
the pressure and the free energy density are connected by the simple relation 

P = • (10) 



On the lattice, the temperature T and the volume V are determined by the 
lattice size x Nr through Eqs. (|l]). The partition function Zl is given by 

Zl{T,V) = y [i5[/]e-'3^(£^^^-^') , (11) 

where U are the link variables and j3A{U\N„,Nr) is the lattice action. 

The free energy density can hardly be determined on the lattice directly by 
a Monte Carlo (which allows to calculate easily expectation values). It can be 
extracted however through the expectation value of the action density St = 
A/ {N^Nr), which is related simply to the derivative of the free energy density 



with respect to the bare coupling P |17|. The free energy density at a given 



temperature T can be thus determined, apart from an additive constant, by an 
integration in /3 up to the value corresponding to that temperature T 



J_ 



rl3 

r4 



/3o JPo 



here the contribution from the action density at zero temperature has been 
subtracted in order to normalize to zero the free energy density at zero temper- 
ature. Actually, (5*0) is calculated on the symmetric lattice iV*, i.e. at a small 
non-zero temperature. Since the low temperature confined phase of SU(N) 
gauge theory is dominated by glueballs which are very massive (mc ~1 GeV), 
the free energy density is expected to drop with the exponential low e"™'^''^. 
This suggests that choosing a small enough /3o fixes the additive integration 
constant to zero (in actual numerical calculations the signal for N^{{Sq) — {St)) 
is zero already at [3 slightly smaller than the critical coupling for the deconfining 
transition). Using Eq. ([l^), the continuous curve for the free energy density as 
a function of the coupling (3 can be determined for a fixed Nt over an arbitrary 
range. The delicate task is to assign a physical temperature T to each value of 
(3 through the relation = Nra{(3). According to the specific lattice action 
adopted, (TN^)^^ may not be small enough in the temperature range of interest 
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to allow the use of the 2-loop asymptotic expression for a{(3). In this case it 
is necessary to improve the determination of a(/3) non-perturbatively, studying 
the /3-dependence of a physical observable^. 

In Refs. [0 numerical simulations were performed using the Wilson action 
with A'^^ = 4, 6, and 8 and N„/N^ > 4. Assuming that N^/N^ was large enough 
to satisfy Eq. @, the curve showing the dependence of p/T'^ on the physical 
temperature T was determined up to temperatures as large as ~ ST^. The non- 
perturbative /3-function was obtained imposing the independence of the critical 
temperature Tc on the critical couplings for Nr = 3,4,6,8 and 12 and using the 
ansatz qAl — f{f3)\2~ioop ■ A(/3), where A(/3) is a smooth function satisfying 
A(/3 — » oo) = 1. It was checked that such a non-perturbative /3-function is in 
agreement with a MCRG analysis of ratios of Wilson loops ||l^ for /3 larger than 
6.0, thus suggesting this value of /3 as the onset of universality. According to the 
critical couplings Pc{Nt) measured in that paper with the Wilson action, only 
for Nt = 8 the curve p/T** for T/Tc > 1 is completely in the universal region of 
the /3-function. This is not true for A^^ = 4, 6 at small T/Tc. 

In Ref. jig] ] a similar program has been carried out using the 2x2 Symanzik- 
improved action on a 24^^ x 4 lattice. The physical scale has been set through 
measurements of the string tension on the 24"^ lattice and of the critical coupling 
on the 24^ x 4 lattice. Interestingly enough, the improvement towards the 
continuum limit occurs even at temperatures slightly above Tc. 

In Fig. ^, the present situation for lattice determinations of the free energy 
density in finite temperature SU(3) gauge theory is summarized^. 



4 Monte Carlo results 



In this Section, the free energy density of finite temperature SU(3) gauge theory 
is determined at temperatures T/Tc = 4/3,3/2 and 2 on lattices with Nr = 
2,3 and No- = ANr, adopting the parametrized FP actions of two different 
RG transformations. The first one is the type I RG transformation defined in 
Ref. ; the second one is the new proposal - type III RG transformation. 
Both actions are parametrized in terms of powers of the real part of the trace of 
the plaquette and of the twisted perimeter-six loop (i.e the loop defined by the 

•^It should be checked, however, that /3 is large enough to prevent any dependence on the 
choice of the observable, i.e. the non-univorsality of the /3-function. 

*I arnjndebted to B^Beinlich, G. Boyd and E. Laermann for giving me tables of the results 
in Ref. M and Refs. M] necessary to build the figure. 
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path xyz-x-y-z, for example). Given the following form for the parametrization 
^(^) = ^ E c,:(C)[iV-ReTr([/c)r , (13) 

C, i>l 

where the sum is over all the closed paths which define a plaquette or a twisted 
loop and Uc = ^cUf_i{n), the couplings Ci(C), i = 1, ... ,4, are determined in 
order to approximate the FP action on a large set of configurations (for details, 
see Refs. |l5[ |6|). They are summarized in Table 1. 



Table 1: Couplings of the parametrized FP action (|13|) for the RG 
transformations of type I ||l5| and type III |l6t]. 







Cl 


C2 


C3 


C4 


type I 


plaquette 


0.523 


0.0021 


0.0053 


.0167 




twisted 


0.0597 


0.0054 


0.0051 


-0.0006 


type III 


plaquette 


0.4822 


0.2288 


-0.1248 


0.0228 




twisted 


0.0647 


-0.0224 


0.0030 


0.0035 



According to Eq. (|l2|), the lattice determination of the free energy density 
at a temperature T calls first for the calculation of expectation values of the 
action density on lattices and x Nr for a sufficiently large set of /3's and 
then for an integration up to the value of f3 corresponding to the temperature 
T at a fixed Nr. 



The usual procedure to determine the correspondence between the coupling 
p and the temperature T at a fixed N-r needs a non-perturbative determination 
of a(/3) or, equivalently, of the /3-function. This procedure usually involves an 
interpolation of numerical data (see the previous Section). Since the central 
point of this paper is to single out lattice cut-off effects of FP actions, an alter- 
native way which introduces no additional artificial effects has been exploited. 
The idea is that the critical coupling PdN^) on a lattice with N!^ sites in the 
time direction and infinite spatial volume corresponds to the physical temper- 
ature T = {Nr/Nr)Tc on a lattice with N^- sites in the time direction. Indeed, 
starting from the definition of the critical coupling /3c (-/V^) 

= iv>(/3c(iv;)) ' ^^^^ 

it is easy to see that it corresponds to the temperature 

1 N'^ 

^ ^ NMMK)) ^ n;^" ^^^^ 

on a lattice with temporal extent Nr- The advantage of this procedure is that 
critical couplings can be determined with high precision. The drawback is, of 
course, that only discrete sets of temperatures can be determined. 
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For the type I FP action, the critical couphngs on lattices with Nr = 2, 3, 4, 6 
were already determined in Ref . ; for the type III FP action the same work 
has also been done In Table 2, results for both actions have been summa- 
rized. 



Table 2: Critical couplings at finite volume and extrapolated to infinite 
volume for the FP actions type I and III parametrized as in Table 1. 





volume 


Nr^2 


Nr = 3 




iV^ = 6 


type I 


43 


2.89(1) 










63 


2.945(15) 


3.320(15) 








83 


2.96(1) 


3.34(1) 


3.50(1) 






10^ 


2.960(5) 


3.34(1) 


3.51(1) 


3.78(2) 




123 


2.975(5) 


3.343(5) 


3.520(5) 


3.775(5) 




153 








3.777(5) 




183 








3.790(5) 




infinite 


2.976(5) 


3.345(5) 


3.520(5) 


3.790(6) 


type III 


43 


3.361(5) 










53 


3.378(3) 










63 


3.385(9) 


3.568(4) 








83 


3.395(3) 




3.678(3) 






93 




3.581(4) 








103 


3.399(5) 




3.686(3) 


3.91(4) 




123 




3.587(5) 


3.691(5) 


3.87(1) 




143 






3.691(5) 


3.882(7) 




163 








3.882(8) 




infinite 


3.400(3) 


3.588(4) 


3.695(4) 


3.885(13) 



According to Eq. ( [l5|) and making use of infinite volume extrapolations in 
Table 2, the f3 couplings corresponding to temperatures the T/Tc ~ 3/2, 2 and 3 
on a lattice with Nt = 2 and those corresponding to temperatures T/Tc = 4/3 
and 2 on a lattice with = 3 can be singled out, for both type I and type III 
FP actions. 

In this paper, simulations have been performed on lattices with — 2 and 
3 and with the ratio Nc/Nr fixed to 4 in order to rule out finite size effect^ 
The simulation algorithm adopted was a mixture of four sweeps of a 20-hit 
Metropolis and one sweep of an over-relaxation consisting in four updates of 
SU(2) subgroups. With the same algorithm, the Wilson action would take ^ 7 
times less for a sweep [ p^ . 

^It is known from measurements of the energy density in SU(2) |Q that thermodynamic 
observables depend on the physical volume through V^/^T = Nc^/Nt-. 
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The free energy density was determined according to Eq. (|l^) at temper- 
atures T/Tc = 3/2 and 2 on the lattice with Nt ~ 2 and at temperatures 
T /Tc = 4/3 and 2 on the lattice with Nr = 3. This allows to compare the re- 
sults at T/Tc = 2 obtained on lattices with different Nr and to check the amount 
of cut-off effects. In Figs. |, | and |, the quantity AS* = iV^((S'o> - (St)) for a 
large range of P is plotted for the Wilson action with N^. — 2 and for the type 
III FP action with N^. = 2 and 3, respectively. In the case of the type I FP 
action the qualitative behavior is the same as for the type III FP action. In all 
three cases the action density shows a rapid variation at the critical coupling, 
which is smoother on the smaller lattices. The time history of the runs for (St) 
in the transition region clearly shows flips between two phases, as expected for 
a first order phase transition. 

For the calculation of the free energy density, AS' has to be integrated with 
respect to f3. For this purpose, a spline fit interpolation of the data for AS* was 
first performed. As an estimate of the uncertainty in the integration coming from 
the interpolation procedure, a straight line interpolation was also considered. 
Although careful checks have been done in the time histories of the runs for St 
in the transition region, the possibility of a (small) systematic error induced by 
hysteresis effects cannot be excluded. In Table 3 all the results for —f/T'^ found 
in this work are summarized. 



Table 3: Values of //T"* at physical temperatures T/T^ = 4/3, 3/2, 2 on 
lattices with Nr = 2 and 3 for the type I and III FP actions and with Nr = 2 

for the Wilson action. 





lattice 


T/T, = 4/3 


T/T, = 3/2 


T/T, = 2 


type I 


8^ X 2 




1.2226(6) 


1.528(4) 




12-3 X 3 


0.769(10) 




1.350(10) 


type III 


8^ X 2 




1.0805(23) 


1.3279(24) 




123 X 3 


0.733(11) 




1.221(10)) 


Wilson 


83 X 2 




2.2007(8) 


2.5652(7) 



In Fig. |6| the results for the type III FP action are compared with pre- 
vious determinations with the Wilson action and the Symanzik-improved 
action , together with the continuum limit extrapolation Q . The horizontal 
error bars represent the uncertainty in the determinations of the physical tem- 
perature T/Tc coming from the uncertainty in the measurements of the critical 
couplings. They have been estimated imposing the independence of the criti- 
cal temperature T^ on the critical couplings for the various A^^-'s and using the 
ansatz aA^ = f{l3)\2-ioop • A(/3), where A(/3) is a smooth function satisfying 
X{P^oo) = 1. 

At A^^ = 2 there are results for -f/T'^ at T/T^ = 3/2 and T/T^ = 2: while the 
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Wilson action shows 0(100%) cut-off effects, tlie type III FP action gives values 
differing from the continuum values estimated in by ~ 24% at T/Tc = 3/2 
and by ~ 9% at T/Tc = 2. 

At Nr = 3, the result for -//T"* with the type III FP action at T/Tc = 4/3 
differs from the continuum by ^ 10% and at T/Tc ~ 2 is consistent with the 
continuum value. 

Although the global scenario exhibits an evident improvement with respect 
to the approach with the standard actions, it is worth to make some comments 
on the surviving cut-off effects in the results of the FP actions. The qualitative 
understanding of such cut-off effects is straightforward. Recalling the classifi- 
cation at the end of Section 2, the effect mentioned in i. should explain the 
deviations between Nr = 2 and Nr = 3 results; the effect mentioned in ii. can 
be responsible of the cut-off effects at small (3 - and so at small T - for each Nr; 
the effect mentioned in iii. can overlap to the previous one at small /3, where 
the parametrization of the FP action is more likely to fail to reproduce the true 
FP action. It would be interesting to evaluate quantitatively the contribution 
of the different sources of cut-off effects: if one of them is largely dominating, 
it can be removed by an ad hoc strategy. 
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Figure 2: The free energy density / versus T/Tc on lattices 16^ x 4, 32^ x 6 
and 32'^ x 8 with the Wilson action (curves labeled with 4, 6 and 8) and 
on lattice 24"^ x 4 with the 2x2 Symanzik-improved action (curve labeled 
with 4(2x2))- The lower curve is the extrapolation to Nt- — > oo obtained from 
the curves with Nr = 6 and 8. The horizontal dashed line indicates the ideal 
gas continuum value. 
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Figure 3: AS* = N^{{Sq) — (■S't)) versus (3 for the Wilson action on a lattice 
8"^ X 2 lattice. The vertical lines represent the critical couplings and the related 
errors on lattices S'^ x 2 [Q, 12^ x 3 |^ and oo^ x 4 0, respectively from left 
to right. 
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Figure 4: AS" = Nf{{Sn) - {St)) versus /3 for the type III FP action on a 8^ x 2 
lattice. The vertical lines represent the critical couplings and the related errors 
on lattices 8^ x 2, oo^ x 3 and oo^ x 4, respectively from left to right (see Table 2). 
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Figure 5: AS = N^{{So) - (St)) versus /3 for the type III FP action on a 12^ x 3 
lattice. The vertical lines represent the critical couplings and the related errors 
on lattices 12^ x 3, oo^ x 4 and oo^ x 6, respectively from left to right (see 
Table 2). 
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Figure 6: The same as in Fig. |^ with the addition of the results of this work for 
the Wilson action and the type III FP action (see Table 3). The vertical dotted 
lines corresponding to T/Tc ~ 4/3, 3/2 and 2 are drawn to guide the eyes. 
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